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Abstract. A precessing source frame, constructed using the Newtonian orbital 
angular momentum Ln, can be invoked to model inspiral gravitational waves 
from generic spinning compact binaries. An attractive feature of such a precessing 
convention is its ability to remove all spin precession induced modulations from 
the orbital phase evolution. However, this convention usually employs a post- 
Newtonian (PN) accurate precessional equation, appropriate for the PN accurate 
orbital angular momentum L, to evolve the L^-hased precessing source frame. 
This influenced us to develop inspiral waveforms for spinning compact binaries 
in a precessing convention that explicitly employ L to describe the binary orbits. 
Our approach introduces certain additional 3PN order terms in the evolution 
equations for the orbital phase and frequency with respect to the usual Ln- 
based implementation of the precessing convention. We examine the practical 
implications of these additional terms by computing the match between inspiral 
waveforms that employ L and Z/^-based precessing conventions. The match 
estimates are found to be smaller than the optimal value, namely 0.97, for a 
non-negligible fraction of unequal mass spinning compact binaries. 
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1. Introduction 

Inspiralling compact binaries containing spinning neutron stars and (or) black 
holes (BHs) are key sources for the network of second generation interferometric 
gravitational wave (GW) detectors [T]. These instruments include the two advanced 
LIGO (aLIGO) observatories [5], the advanced Virgo the KAGRA [3], the GEO- 
HF [5] and the planned LIGO-India [6]. In contrast, massive spinning BH binaries 
are one of the most exciting sources for the space-based GW observatory like the 
planned eLISA [7j. GWs from such inspiralling compact binaries, whose components 
are specified by their masses and spins, can be accurately modeled using perturbative 
approaches to tackle the underlying Einstein field equations [5]. Therefore, the 
optimal detection technique of matched filtering can be employed to detect and 
characterize inspiral GWs from such binaries [anQ]. This involves cross correlating 
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the interferometric output data with a bank of templates that theoretically model 
inspiral GWs from spinning compact binaries. A successful detection demands that 
at least one template should remain in phase as much as possible with the buried weak 
GW signals in the frequency windows of various GW observatories m- 

During the GW emission induced inspiral, dynamics and associated GWs 
from compact binaries can be accurately described using the post-Newtonian (PN) 
approximation to general relativity mm- The PN description provides interesting 
quantities, required for various template constructions, as an asymptotic series in 
terms of certain dimensionless parameter. For binaries in quasi-circular orbits, it is 
usual to use x = (Gmw/c^)^/^ as the PN expansion parameter while constructing 
inspiral templates, where m and uj stand for the total mass and orbital (angular) 
frequency of the binary [n m. Gurrently, GW frequency and associated phase 
evolution, crucial inputs to construct various inspiral template families, are known 
to 3.5PN order for non-spinning compact binaries [m US]- In other words, PN 
corrections to the above two quantities are computed to the order beyond the 
leading quadrupolar (Newtonian) order for such binaries. Moreover, the amplitudes of 
the two GW polarization states, hx and /i+, for non-spinning binaries are available to 
3PN order m- Very recently, detailed computations led to the determination of the 
dynamics of such binaries to the 4PN order [18] . Binaries that contain compact objects 
with intrinsic rotations, the spin effects enter the dynamics and GW emission via spin- 
orbit and spin-spin interactions naiio]. In binaries containing maximally spinning 
BHs, the spin-orbit coupling (linear in the spins) first appears at the 1.5PN order, 
while the spin-spin interaction (which is quadratic in spins) first occurs at the 2PN 
order |2T]. Additionally, Si, S 2 and L, the two spin and orbital angular momenta, 
for generic spinning compact binaries precess around the total angular momentum 
J = L Si S 2 due to spin-orbit and spin-spin interactions. This forces substantial 
modulations of the emitted GWs from inspiralling generic spinning compact binaries 
compared to their non-spinning counterparts [211 [22]. Therefore, it is important to 
incorporate various effects due to the intrinsic rotations of compact objects while 
constructing inspiral GW templates for spinning compact binaries. At present, GW 
frequency evolution and amplitudes oi hx and for BH binaries having maximally 
spinning components are fully determined to 2.5PN and 2PN orders, respectively, while 
incorporating all the relevant spin induced effects [lllESlEllEsJEelET]. Moreover, 
the on-going detailed computations are providing various higher PN order spin-orbit 
and spin-spin contributions to the dynamics of spinning compact binaries in general 
orbits and to the orbital frequency evolution for quasi-circular inspiral. At present, 
the spin-orbit contributions to binary dynamics and GW frequency evolution are 
available up to the next-to-next-to-leading order (2PN order) beyond the leading order 
[28l [29] while adapting the MPM (Multipolar post-Minkowskian) approach [30]. In 
contrast, the higher order spin-spin contributions to the binary dynamics are usually 
tackled in the Arnowitt-Deser-Misner canonical formalism m and in the Effective 
Field Theory formalism [32] [33] (Note that the spin-orbit effects in the Effective Eield 
Theory formalism are computed, for example, in [34]). These approaches provide 
various spin(l)-spin(2) and spin-squared contributions to the orbital dynamics |35j . 
Moreover, various source multipole moments needed to obtain the spin contributions 
to GW luminosity at the 3PN order and GW polarization states to the 2.5PN order 
were computed in [55] . 

The construction of time-domain hx and /i+ associated with inspiralling generic 
spinning compact binaries requires us to numerically solve a set of PN accurate 
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differential equations for Si, 82 , L^,x and the orbital phase m, where Xn is the 
Newtonian orbital angular momentum. The numerical integration provides temporal 
evolutions for the orbital phase, the associated angular frequency and the two angles 
that specify the orientation of the orbital plane in an inertial frame associated with 
the direction of J at the initial epoch. These variations are incorporated into the 
PN accurate expressions for hx and /i+ to obtain PN accurate time-domain inspiral 
waveforms from such binaries |21] . In this approach, the differential equation for 
the orbital phase explicitly depends on the precessional motion of the orbital plane 
mull]. Therefore, it is not possible to express the orbital phase as an integral of the 
orbital frequency as usually done in the case of non-spinning compact binaries m- 
A decade ago, Buonanno, Chen and Vallisneri proposed an approach, referred to as 
the precessing convention, that factorizes the generic spinning binary waveform into 
a carrier signal and a modulated amplitude term m- In this approach, the phase of 
the carrier signal (<I>p) essentially coincides with the accumulated orbital phase such 
that <I>p = uj. Moreover, the precessional dynamics of the orbital plane only influences 
the modulated amplitude part of inspiral waveform even for generic compact binaries. 
Therefore, the precessing convention disentangles the precessional effects from its non- 
precessional counterparts while modeling both the amplitude and the phase evolutions 
of inspiral GWs from such astrophysical systems. This convention was employed to 
model inspiral GW signals from compact binaries containing misaligned single-spin 
and to probe its data analysis benefits [381139] . Very recently, inspiral-merger-ringdown 
waveforms for generic spinning BH binaries, invoking the effective-one-body approach 
[40] . also adapted the precessing convention to model GWs from the inspiral part [4T] . 
We note that this convention requires a precessing source frame which is usually based 
on the Newtonian orbital angular momentum Tn- However, [41] employed both L 
and Ljq to model GWs during the late part of the binary inspiral just prior to the 
plunge. 

In this paper, we develop a prescription to compute PN accurate inspiral 
waveforms for generic spinning compact binaries while using the PN accurate orbital 
angular momentum L to construct the precessing source frame. This is influenced 
by the usual practice of employing precessional equation appropriate for L to evolve 
Tn and the associated precessing source frame while constructing inspiral waveforms 
via the precessing convention of [37] • We show that the use of such an adiabatic 
approximation, namely employing an orbital averaged differential equation for Tn, 
can lead to PN corrections to = ui. These observations motivated us to provide 
a set of PN accurate equations to obtain temporally evolving quadrupolar order 
inspiral GW polarization states for generic spinning compact binaries in the X-based 
precessing convention. In our approach, the spin precession induced modulations 
enter the differential equation for the orbital phase only at the 3PN order. Moreover, 
the X-based convention requires us to include additional 3PN order terms in the 
differential equation for x, compared to the usual XN-based approach. We explore the 
practical implications of these additional terms with the help of match computations, 
detailed in [ni Eg. The match computations involve two inspiral families where 
one is constructed via our X-based precessing convention and therefore incorporate 
the above mentioned 3PN order terms. The other family is based on the lalsuite 
SpinTaylorT4 code, developed by the LIGO Scientific Gollaboration (LSG) [43], that 
implemented the XN-based precessing convention of [37] . In this code, the precessional 
dynamics is fully 2PN accurate while the differential equation for x incorporates spin- 
orbit contributions to 3.5PN order. For match computations, we employ PN accurate 


On precessing convention 


4 


relation connecting and L to construct two waveform families with physically 
identical initial orbital and spin orientations. We terminate the two inspiral waveform 
families when their respective x parameter reaches 0.1 which roughly corresponds 
to orbital separations ~ 10 Gm/c^. These considerations allow us to attribute the 
reported match estimates to the above mentioned additional 3PN order terms present 
in our differential equations for the orbital phase and frequency. We find that the 
match estimates are less than the optimal 0.97 value for a non-negligible fraction of 
unequal mass spinning compact binaries. It may be recalled that such an optimal 
match value roughly corresponds to a 10% loss in the ideal event rate. 

In what follows, we briefly summarize the usual implementation of the L^- 
based precessing convention and explore the consequence of employing an orbital 
averaged precessional equation for Tn. In section 12.21 we detail the construction 
of quadrupolar order GW polarization states in our T-based precessing convention. 
The match estimates involving these two inspiral families with physically equivalent 
initial configurations and associated discussions are listed in section |3] while section!?] 
provides a brief summary. 

2. Inspiral 'waveforms in and i-based precessing conventions 

We begin by summarizing the usual implementation of the iN-based precessing 
convention and explore the consequence of employing an orbital averaged differential 
equation for Tn- 

S.L The -based precessing convention 

The precessing convention, introduced in [37], aims to remove all the spin precession 
induced modulations from the orbital phase evolution. In this approach, the orbital 
phase ^p{t) is written as an integral of the orbital frequency uj{t), namely $p(t) = 
f uj{t') dt', even for generic spinning compact binaries. This feature is crucial to ensure 
that the inspiral waveform for a precessing spinning compact binary can be written as 
the product of a non-precessing carrier waveform and a modulation term that contains 
all the precessional effects. This is how the approach disentangles the precessional 
effects from their non-precessional counterparts both in the amplitude and phase of 
inspiral waveforms. It should be noted that in the absence of precessing convention, 
the orbital phase of a generic spinning binary is given by J[u}{t') — a'{t') cosi'(t')] dt', 
where i' and o' specify the orientation of Tn in an inertial frame associated with the 
initial direction of J |2I) . 

To obtain inspiral waveforms for spinning compact binaries in their precessing 
convention, m employed certain precessing source frame = 1 ), where I is 

the unit vector along Tn ■ The basis vectors of this triad satisfy the evolution equations 
^ 1.2 = ^ ^ 1,2 nnd e\ = I = rik 1. The angular frequency Vt\ is constructed in 

such a manner that these three basis vectors always form an orthonormal triad. This 
is possible with the following expression for fig, namely fig = fi^ — (fi^ • 1 ) 1 , where 
fifc is the usually employed precessional frequency for 1. The expression for fi^ that 
includes the dominant order spin-orbit and spin-spin contributions can be obtained by 
collecting the terms that multiply Tn (= 1) in equation (9) of [37]. With the help of 
e[ and e^, an orbital phase ^p{t) may be defined such that n = cos e\ -\- sin <I>p 
where n is the unit vector along the binary separation vector r. Additionally, one may 
define a co-moving frame (n, X — lxn,l) such that the time derivative of n is given by 
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n = $pA+r2g xn. It was argued in m that rig should only be proportional to n which 
ensures that n - A = 4)p. This leads to the desirable expression for 4)^, namely $p = w, 
while employing the adiabatic condition for the sequence of circular orbits: h - \ = to. 
This adiabatic condition is equivalent to n ■ h = uP' that provides another version 

of the PN independent relation connecting the linear and orbital angular velocities, 

2 2 2 
V = u) . 

A close inspection reveals that this approach usually employs an orbital averaged 
differential equation for L-^ to evolve the precessing source frame while constructing 
PN accurate inspiral waveforms. It turns out that the differential equation for L-^ 
in such an adiabatic approximation is identical to the evolution equation for the PN 
accurate orbital angular momentum L [3310^. In what follows, we explore the effect 
of such an adiabatic approximation on the equation for fi in the (n, A,Z) co-moving 
frame and on the derivation of <I>p equation. The usually employed expression for 
while considering only the leading order spin-orbit interactions may be written as 


= 7^ a: 
Gm 


^ ^ hqxisi + —X2S2 


( 1 ) 


where q = rfixjrfi^ (mi > m 2 ) is the mass ratio and ^ 1^2 = ??/2 -I- 3 (1 =F ~ 4?7)/4 
while T] = mi m^lrrp is the symmetric mass ratio. The Kerr parameters xi and X 2 of 
the two compact objects of mass mi and m 2 specify their spin angular momenta by 
S \_2 = Gmf 2 X 1.2 where Si and S 2 are the unit vectors along Si and S 2 - It is 

straightforward to show that the above equation is identical to terms that multiply 
Ln on the right hand side of equation (9) in [37j . To explore the implication of using 
the above expression for flk to construct Jig, we revisit the arguments detailed in the 
appendix B of m- These arguments, crucial to obtain 4>p = w, require that fig x n, 
appearing in the equation n = <i>pA-|-fig x n, should be zero. A closer look reveals that 
m did not use the explicit expression for fig to show that fig lies along n. Instead, 
the authors arrived at such a conclusion with the help of the following two steps. First, 
it was noted that I should be oc A (see lines around equations (B3), (B4) and (B5) in 
[37]). This was inferred by employing the definition for Z (Z = n x A), the resulting 
time derivative for Z {I = nx\ + nx\), the adiabatic condition for circular orbits 
(n = w A) and the time derivative for A in the co-moving frame (A = — 4>p n-t-flg x A). 
In the second step, [37] invoked the requirement that Z should also be given by fZg x Z in 
the precessing source frame. With the help of the above two arguments, namely Z oc A 
and Z = fZg X Z, [S^ concluded that fZg should lie along n as \ — lxn. We would like 
to emphasize that |37j never invoked their explicit expression for fZg to show that fZg 
can have components only along n which is essential to obtain the relation $p = w. 
In fact, it is straightforward to show with the help of equation o that fZg can have 
components along A since fZg • A = /Gm) (di qxi^i ■ X + 82 X 2 /qs 2 • A) ^ 0, in 
general. This results in the following 1.5PN accurate expression for fi 

n = 4>pA -F ^ gxi [si xn - {si • Z) A] -f ^ X 2 [^2 x n - (s 2 • Z) A] | . (2) 

Clearly, the 1.5PN order terms that arise from fZg x n in the above equation will not 
be zero for generic spinning compact binaries. 

Interestingly, it is still possible to equate 4>p to w by employing the adiabatic 
condition n • A = w as (fig x n) • A = 0 even in the presence of non-vanishing 1.5PN 
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order fig x n term. However, the equivalent version of the adiabatic condition, namely 


h ■ n = uj^, forces the differential equation for to contain 3PN order corrections 
in addition to the standard w term. In other words, two equivalent versions of the 
same adiabatic condition for circular orbits, namely h ■ X = uj and fi ■ fi = 
provide different evolution equations for $p. This is the unexpected consequence of 
employing processional equation appropriate for L to evolve the TN-based precessing 
source frame. It is not difficult to deduce that this essentially arises from the non¬ 
vanishing I.5PN order ft], x n contributions to n listed in the above equation. The 
following arguments can also be used to clarify why the above two versions of the 
same adiabatic condition for the circular orbits result in two different expressions for 
<I>p. In our opinion, this arises by identifying k x n to he X = I x n, where k is the 
unit vector along L. Strictly speaking, the precessing source frame of m is based on 
L rather than Tn due to the use of that provides the precessional equation for L 
|441145) . This implies that their co-moving triad is rather i-based and the expression 
for fi got components along k x n instead oi X = I x n. Therefore, n • A results in 
PN corrections to as (fc x n) • A is unity only at the leading order (this may be 
deduced from our equation ([ 8 ]) listed below). This leads to a differential equation for 
$p that involves PN corrections to ui. 

The use of precessional equation appropriate for L while implementing the 
precessing convention of m motivated us to develop a fc-based precessing convention 
for constructing inspiral waveforms for spinning compact binaries. This should also 
allow to explore the practical implications of using the adiabatic approximation to 
evolve Tn in the usual implementation of the precessing convention. 

2.2. Inspiral waveforms via an L-based precessing convention 

Influenced by the above arguments and m, we first introduce a fc-based precessing 
source frame: (ei,e2, 63 = fc). The precessional dynamics of ei, 62 and 63 are provided 
by 61^2,3 = X 61^2,3, where He = — (H^-fc) fc and is the precessional frequency 

of fc. It should be obvious that 63 = He x 63 is identical to 63 = H^ x 63 as 63 = fc. 
It is possible to construct a fc-based co-moving triad — k x n, fc) and introduce 
an orbital phase $ such that 


(3a) 

(3b) 


n = cos $ ei -I- sin $ 62 , 

^ = — sin $ ei -I- cos $ 62 . 


It is fairly straightforward to obtain following expressions for the time derivatives of 
n and ^ 


(4a) 

(4b) 


n = $^-|-HeXn, 

^ = —$ n -I- He X ^ . 


We are now in a position to obtain the differential equation for $. This is derived 
with the help of the frame independent adiabatic condition for circular orbits, namely 
h-n = u!^. Employing the above expression for h in such an adiabatic condition leads 
to 



(5) 


where Hej = He • ^ is given by 



(6) 
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This results in the following 3PN accurate differential equation for $ 


$ = 


p3/2 


Gm 


1 - 


r ^2 

^19X1 («1 -i) + — X2(S2 -C) 


2 L 


(7) 


These additional terms appear at the 3PN order as we employ equation ([IJ for $7^ that 
only incorporates the leading order spin-orbit interactions appearing at 1.5PN order. 
It is not difficult to deduce that the inclusion of 2PN order spin-spin interaction terms 
in Ofc can lead to certain 4PN order contributions to the $ equation. Additionally, 
the presence of the above 3PN contributions to $ equation may be attributed to the 
fact that the orbital velocity v can have non-vanishing PN order components along L 
while considering generic spinning compact binaries [44) . 

The use of L to describe binary orbits also modifies the evolution equation for uj 
(or x). This is because the spin-orbit interactions are usually incorporated in terms 
of Si • i and S 2 ■ I, in the literature [ 211137 ]. These terms require modifications due to 
the following 1.5PN order relation connecting I and k 

l^k + ^ [xi («i ■n)+X2 (S2 ■ w)] n 

+ [2 Xi (si • 0 + ?? Xi (si • 0 

+ 2 A 2 X2 {S2 ■ $,) + 11 x 2 (S2 • I)] ^I , (8) 

where Xi = mijm and X 2 = m^lm. The above relation can easily be extracted, 
for example, from equations (6.10) and (7.10) in [22]- This relation leads to certain 
additional 3PN order contributions to x while describing the binary orbits with k. 
These additional terms are, for example, with respect to equation (3.16) of [ 25 ] that 
provides PN accurate expression for dx/dt while invoking I to describe binary orbits. 
Our additional contributions to x appear at 3PN order as the dominant spin-orbit 
interactions, in terms of Si • Z and S 2 ■ Z, contribute to the x evolution equation at 
1.5PN order. The differential equation for x in our approach may be written as 


dx dx 

— = —{Equation (3.16) in [21]; Z —>■ fc) 

(JjL (JjL 


64 

5 Gm 


rjx <x 


47 , 

-y ^1 Xi 


\d(xi (si ■n)+X2 {S2 ■ n)'j (si • n) 


+ (2 X^ xi («i ■ I) + ?? Xi (si ■ 0 + 2 X| X2 (S2 ■$)+dX2 (S2 ■ D) (si ■ ^) 

- y xi X2 V (xi (si ■ ’^) + X2 (S2 ■ n)) (^2 • n) 

+ (2 xi Xl («1 ■ I) + ?? Xl (Sl ■ 0 + 2 xi X2 {S2 ■O+VX2 {S2 ■ D) {s2 ■ 

- y (^1 - X2) X2 X2 (xi (si ■n)+X2 (s2 ■ «)) (s2 • n) 

+ (2 xi Xl (si ■ 0 + ?? Xl (si • 0 + 2 xi X 2 {S 2 ■O+VX2 {S 2 ■ 0) (S 2 ■ $,) 
+ y (^1 - X2) Xl Xl (xi (si ■ »^) + X2 (S2 ■ n)) (si ■ n) 

+ (2 Xi Xl (Sl • 0 + ^ Xl (si • 0 + 2 xi X2 {S2 ■$.)+VX2 {S2 ■ 0) (si • 0 


(9) 
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where the first term is adapted from equation (3.16) in by replacing its I vectors by 
our k vectors. Clearly, these additional spin-squared terms contribute to dx/dt at 3PN 
order. Note that the use of k for I in the leading order spin-spin contributions to dx/dt, 
as listed in equation (1) of |3Z]: results in additional 3.5PN order spin-cubed terms. 
Such contributions to our dx/dt equation are neglected as the lalsuite SpinTaylorT4 
code that implement the iN-based precessing convention does not include any spin- 
cubed terms in its differential equation for x. 

In what follows, we model inspiral GWs from spinning compact binaries in our k- 
based precessing convention. We begin by displaying the following quadrupolar order 
expressions for the two GW polarization states: 

hxlqit) =2 (2 no, Uy ), ( 10 a) 

( 10 b) 

where ^x,y and are the x and y components of ^ and n in an inertial frame 
associated with TV, the unit vector that points from the source to the detector, while 
R' is the distance to the binary. These x and y components can be expressed in terms 
of the Cartesian components of ei and 62 via equations (jS]). We note that the above 
expressions for the quadrupolar order and Ti+ are written in the so-called frame¬ 
less convention m- These expressions emerge from the following standard definitions 
for the quadrupolar order GW polarization states: 

hxlqit) = +q^p’)hJ^\Q, ( 11 a) 

h+\Q{t) = ^ (pV - q"q^)hJ^\Q , ( 11 b) 

where p and q are the two polarization vectors forming, along with TV, an orthonormal 
right-handed triad. To obtain equation we used the following expression for the 
quadrupolar order transverse-traceless part of the far-zone field 



AGmrix 
c2 R' 


(re 


n* n^). 


( 12 ) 


In the frame-less convention, we let the components of p and q in the TV-based inertial 
frame to be p = (1,0,0) and q = (0,1,0). Clearly, we need to specify how the 
Cartesian components of ^ and n vary in time to obtain temporally evolving GW 
polarization states for inspiralling generic spinning compact binaries. Therefore, we 
require to solve numerically the differential equations for $,a:, ei and 62 to obtain 
hxlqit) and Ti+| q(T). We use equation 0 for $ while the differential equation for x 
is given by equation ([ 2 ]) that contains all the non-spinning contributions accurate up 
to 3.5PN order and the usual spin contributions that are fully 2PN accurate. These 
contributions are provided, for example, by equation (3.16) in and are listed in 
the [Appendix A[ where we have replaced I hy k. These specific PN order choices are 
influenced by the fact that the lalsuite SpinTaylorTd code also employs an evolution 
equation for w that incorporates such PN contributions. We note that this routine 
implements the T-based precessing convention of m to construct inspiral templates 
to search for GWs from generic spinning binaries. The differential equations for ei 
and 62 in our approach are given by 

e'l == r2e X ei = {flk - (Ofe • k)) x ei, 

62 = r2e X 62 = {flk “ {^k ' k)) X 62 , 


(I3a) 

(13b) 
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where we use the following expression for that incorporates the leading order 
spin-orbit and spin-spin interactions 


= ^ 2 ;^ \hqXiSi + ^X2S2- X 2 [{k ■ Si) S 2 + {k ■ 82 )^ 1 ] 


(14) 


It is fairly straightforward to verify that this expression is identical to the coefficient 
of Ln appearing on the right hand side of equation (9) in [37] . The above equations 
imply that we also need to invoke the precessional equations for Si, S 2 and k (or 63 ) 
to tackle numerically the dynamics of such binaries. The three coupled equations for 
Si, S 2 and k that include the leading order spin-orbit and spin-spin interactions read 


Si = -I (5i (fc X si) -b ^ a; 

Gm ' ^ 


1/2 


Xix 2 (S 2 X Si) -3X^X2 {k ■ S 2 ) {k x Si) 


(15a) 




ikxS2) + ^X^/^ 


Xi X1 (si X S 2 ) - 3 X^ Xi (fe ■ Si) (fc X S 2 ) 


(15b) 


k = x^ < Si q xi (si X fc) H-X2 (S2 X k) 


Gm 
3 


x^^'^VXi X2 


(k ■ Si) (S 2 X fc) -b (fc • S 2 ) (si X fc) 


(15c) 


It is not very difficult to verify that the above equations for Si and S 2 are identical 
to equations ( 2 ) and (3) in [57] while the differential equation for fc, as expected, 
arises from the usual conservation of total angular momentum J. This conservation 
implies that Lk = —Si s'l — S 2 S 2 . We would like to state again that the equations 
([7]), (lA.ll) . (IA.2P and ([9]) provide the differential equations for $ and x in the 
present implementation of fc-based precessing convention. Strictly speaking, the use of 
equation (IT4ll for flfc requires us to include the additional 4PN and 3.5PN contribution 
to d^/dt and dx/dt, respectively. However, we do not incorporate such spin-quartic 
and spin-cubic terms in our present work. 

In practice, we numerically solve simultaneously the differential equations for 
ei, fc, Si, S 2 i d> and x to obtain temporally evolving Cartesian components of n and 
The resulting variations in these Cartesian components are imposed on the expressions 
for fcx,-i-lQ(^)i given by equations (fTUl) . This leads to inspiral waveforms for generic 
spinning compact binaries in our fc-based precessing convention. Note that we do not 
solve the differential equation for 62 . This is because the temporal evolution of 62 
can be estimated using the relation e 2 (t) = k{t) x ei{t). This implies that we solve 
12 differential equations for the Cartesian components of ei, fc, si and S 2 along with 
the differential equations for $ and x to track the time evolution for the Cartesian 
components of n and The required initial values for the Cartesian components of 
ei, fc, Si and S 2 are given by freely choosing the following five angles: 9iq, c/iiq, 6 * 20 , 
4 >2 o and lq. These angles specify the above four unit vectors in the TV-based inertial 
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frame at the initial epoch such that 

51 = (sin 010 cos(/)io, sin 010 sin(^io,cos0io) , (16a) 

5 2 = (sin 020 cos()) 2 o,sin 02 o sin(/) 2 o,cos 02 o) , (16b) 

fc = (sinio,0,cos6o), (16c) 

ei = (cos tO) 0, — sinio) • (16d) 


This choice is also influenced by the lalsuite SpinTaylorT4 code of LSC. Additionally, 
we let the initial x value to be xo = (Gmcuo/c^)'^^^ where ujo = IOtt Hz and the initial 
phase $0 to be zero. 

We move on to explain how to specify initial conditions, physically identical to 
equations (HSI), while constructing inspiral waveforms based on the Z-based precessing 
convention. Clearly, the initial orientations of two spin vectors in the inertial N frame 
should be identical in the two approaches. However, the orientation of I in such an 
inertial frame is different from that of k. We compute TV • I using equation ([5]) that 
provide 1.5PN accurate relation connecting I and k. This PN relation makes the value 
of TV • Z at Xo to depend on lq, m, t], Xi, X 2 ) ‘f’o and the four angles that specify Si and 
S 2 in the inertial TV frame. We observe that [TV • Z(xo) — TV • fc(xo)] is maximum 
for equal mass maximally spinning compact binaries and the difference is usually less 
than 0.1%. The value of TV • l{xo) specifies the initial orientation of Z in the inertial 
TV frame as we usually let its azimuthal angle to be zero along with $p(xo) (we have 
verified that the changes in these angles play no role in our match computations). 
The difference in TV • 1{xq) and TV • k{xo) values leads to slightly different values 
for fc • Si and Z • Si at the initial epoch. Therefore, the orbital frequency and phase 
evolutions are slightly different in the two approaches even in the absence of our 3PN 
order additional terms. We observe that the differences in these two dot products, 
namely k- Si and Z • Si, have their maximum values for equal mass maximally spinning 
compact binaries. In the next section, we purse the match computations to probe 
the implications of additional 3PN order terms present in the frequency and phase 
evolution equations while constructing our fc-based inspiral waveforms. 

3. Match computations involving the above two families of inspiral 
waveforms 

We employ the match, detailed in [TTJ |42] , to compare inspiral waveforms constructed 
via the above described Z and fc-based precessing conventions. Our comparison is 
influenced (and justified) by the fact that the precessing source frames of the above 
two conventions are functionally identical. This should be evident from the use of the 
same processional frequency, appropriate for fc, to obtain PN accurate expressions for 
both the Z-based Hg and fc-based Hg- Therefore, the match estimates probe influences 
of the additional 3PN order terms present in the differential equations for $ and x 
in our approach. Additionally, we have verified that these 3PN order terms are not 
present in the usual implementation of the precessing convention as provided by the 
LALSUITE SpinTaylorT4 code. 

Our match M{hi, hk) computations involve hi and h^, the two families of inspiral 
waveforms arising from the Z and fc-based precessing conventions. The hi inspiral 
waveform families are adapted from the lalsuite SpinTaylorT4 code of LSC while hk 
families, as expected, arise from our approach (equations (flUl) '). It should be noted that 
we employ the quadrupolar (Newtonian) order expressions for hx,+ while computing hi 
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and hk in the present analysis. However, the lalsuite SpinTaylorT4 code can provide 
waveforms that include 1.5PN order corrections to their amplitudes. We would like to 
stress that the two families involved in our M{hi, hf^) computations are characterized 
by identical values of to, 77 , xi and X 2 - Additionally, the initial orientations of the two 
spins in the Af-based inertial frame were also chosen to be identical. The computation 
oi N ■ I from N ■ k with the help of equation ([8]) ensures that I and k orientations 
at the initial epoch are physically equivalent. Therefore, our Ai{hi, hk) computations 
indeed compare two waveform families with physically equivalent orbital and spin 
configurations at the initial epoch. To obtain a specific M.{hi,hk) estimate, we first 
compute an overlap function between the relevant hi and hk inspiral waveforms: 


0{hi,hk) =< hi,hk >= 


{hi\hk) 

^/{hl\hl) {hk\hk) ’ 


(17) 


where hi and hk stand for the normalized hi{t) and hk{t) waveforms, respectively. 
The angular bracket between hi and hk defines certain noise weighted inner product, 
namely 


{hi\hk) = 4Re 



hnf)hk{f) 

SM) 


(18) 


In the above equation, hi{f) and hk{f) stand for the Fourier transforms of hi(t) 
and hk{t), while S'h(/) denotes the one-sided power spectral density. We have 
invoked the zero-detuned, high power sensitivity curve of aLIGO [35] in om M{hi, hk) 
computations. The upper cut-off frequency /cut is chosen to be c^/(G to tt 10^/^) while 
the lower cut-off frequency /low, associated with the GW detector, equals 10 Hz. The 
match A4(hi, hk) is computed by maximizing the 0 {hi, hk) over two extrinsic variables, 
namely the time of arrival <0 and the overall phase (po of GW at time to jUj- This 
leads to 

Ai = max O{hi, hk) ■ (19) 

to,4>o 


The maximizations of over to and 0o performed by following m- We perform the 
maximization over to via the FFT algorithm while the maximization over pg requires 
us to deploy two orthogonal templates. Let us emphasize that we terminate hi and 
hk inspiral waveform families when their respective x parameters reach 0.1 which 
roughly corresponds to orbital separations ^ 10 Gto/c^. This choice arguably ensures 
the validity of PN approximation to describe the temporal evolutions of the above two 
families in our [/iow-/cut] frequency window. Therefore, it is reasonable to associate 
the departure of M. estimates from unity to the additional 3PN order contributions to 
the differential equations for the orbital phase and the associated angular frequency. 

We move on to list results of some of our M. computations. It is not very 
difficult to realize that extensive M. computations that deal with all aLIGO relevant 
spin and binary configurations will be rather difficult to achieve. Therefore, we 
restrict our attention to a selected number of binaries to compare the hi{t) and 
hk{t) inspiral waveform families. In our match computations, we mainly consider 
binaries having total mass to > 30Mq due to the following two reasons. First, it is 
comparatively expensive (computationally) to generate lengthy inspiral waveforms for 
low mass binaries in the aLIGO frequency window. Additionally, we are interested 
in to explore the dependence of our Ai estimates on the mass ratio q in the [1 — 10] 
range. Clearly, low to binaries can lead to secondary BHs having masses lower than 
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the usual neutron star masses for high q cases. In figure [TJ we consider maximally 
spinning BH binaries, characterized by the total mass m = SOMq, while varying the 
mass ratio q from unity to 10. The left and right panel plots are for configurations 
having initial dominant spin-orbit misalignments 0i(a:o) (cos“^(fc • Si)) given by 30° 
and 60°, respectively. Additionally, we let the initial orbital plane orientation in 
the Af-based inertial frame to take two values leading to edge-on {lq = 90°) and 
face-on (to = 0°) binary orientations. These binary configurations, characterized by 
two different spin-orbit misalignments and orbital plane orientations in the TV-based 
inertial frame, are obtained by choosing appropriately different values for the more 
massive BH initial spin orientation, namely 0i(xo) in equation (I16al) . For example, we 
choose 010 to be 30° and 60° to get 9i{xo) equals 30° and 60°, respectively, for face-on 
binaries. However, in the case of edge-on binaries, we require to choose 0io to be 60° 
and 30°, respectively. All other initial angular variables, appearing in equations (I16ap 
and (I16bl) . were chosen to be ^lo = 0°, 02O = 20°, (j) 2 o = 90° (we have verified that the 
Ai estimates are rather insensitive to the choice of these initial angular variables). Let 
us note that I orientations (from TV) for these configurations will be slightly different 
from 0° or 90° due to the 1.5PN accurate relation between I and k. 






Figure 1. Plots for the accumulated orbital phase (A$) and the associated match 
(Ad) estimates as functions of the mass ratio q for maximally spinning m = SOMq 
BH binaries inspiralling in the [/iow"/cut] frequency interval. The left and right 
panel plots are for binaries having dominant spin-orbit misalignments 6 i = 30° 
and 60°, respectively, at xq. The solid line plots in each panel correspond to 
face-on binaries whereas the dotted line curves are for edge-on binaries. We find a 
gradual degradation of the Ad estimates for higher q values for 6 i = 30° binaries. 
However, this trend is not visible for 0i = 60° binaries. These configurations 
with different spin-orbit misalignments and binary orientations are obtained by 
setting different values for ^lo and lq. Face-on binaries with 6 i{xq) = 30° and 60° 
have (^ 10 ,^ 0 ) ^ (30°,0°) and (60°,0°) while the edge-on binaries have (^ 10 ,^ 0 ) 
as (60°, 90°) and (30°, 90°), respectively. The initial values of the other angular 
parameters remain the same, namely 0io = 0°, 620 = 20°, 020 = 90°. The 
PN accurate relation between h and I also results in slightly different dominant 
spin-orbit misalignments in the two approaches. 

The upper and lower row plots in figure [1] are for A$, the accumulated orbital 
phase differences in the [/iow-/cut] frequency interval while dealing with the hi and 
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hk inspiral waveforms and the associated Ai{hi,hi^) estimates, respectively. We find 
that the variations in M estimates are quite independent of the initial orbital plane 
orientations. However, variation in match estimates do depend on the initial dominant 
spin-orbit misalignment and the mass ratio q. The left panel plots show gradual 
decrease in M values as we increase the q value and this variation is reflected in the 
gradual increase of A<i>. Incidentally, this pattern is also observed for configurations 
having somewhat smaller initial dominant spin-orbit misalignments. However, the M 
estimates are close to unity for tiny 0 i{xo) values and this is expected as precessional 
effects are minimal for such binaries. Therefore, the effect of the above discussed 
additional 3PN order terms are more pronounced for high q compact binaries having 
moderate dominant spin-orbit misalignments. This picture is modified for binaries 
having substantial dominant spin-orbit misalignments as evident from the M plots 
in the right panel of figure [I] For such binaries, the M estimates dip to a minimum 
and recover as we vary q from 1 to 10 for both edge-on and face-on orbital plane 
orientations. In contrast, we observe a gradual increase in A^{q). The monotonic 
increment in A<I>(( 7 ) plots is essentially due to the presence of Xi — X 2 terms in the 
additional 3PN order contributions to dx/dt, given by equation (jH). This is because 
Xi — X 2 terms are absent for equal mass binaries which leads to smaller A$ estimates 
for smaller q value binaries. However, it is not possible to explain the observed M{q) 
variations purely in terms the displayed A^(q) values, especially when precessional 
effects are substantial as in the 9i{xo) = 60° cases. We observe that the initial values 
of angles like cos“^(fc • ji) and cos“^(Ar ■ j) also influence the precessional modulations 
in the waveforms, where j is the unit vector along J. Therefore, we examined how 
the initial values of these angles vary as functions of q. We find monotonic increase 
(decrease) in the initial values of cos~^ {k ■ j) {cos~^{N -j)) when q value is varied from 
I to 10. Therefore, we speculate that the combined effect of such variations and the 
non-negligible A$(g) values may provide a possible explanation for the dip in match 
estimates around q=5. This is because initial values of the above mentioned angles do 
define the way various dot products, involving n, 4, Si and S 2 , vary in time. It should 
be noted that these dot products are present in the additional 3PN order terms in the 
dx/dt expression given by equation (jH]). 

This plausible explanation is tested in figure [2l where we plot A$(g) and M[q) 
estimates for binary configurations having three different m values. It is not difficult to 
infer that the above listed arguments should also hold for such compact binaries. The 
initial dominant spin-orbit misalignment is again chosen to be 60° while considering 
only edge-on configurations. All other spin and orbital orientations at the initial epoch 
are identical to the cases displayed in the right panel plots of figure [TJ Clearly, the 
plots in figure [2] are qualitatively similar to those in the right panel plots of figure [I] 
However, the dip in M. estimates shifts to higher q values for higher m compact 
binaries. We conclude from the above two figures that GW data analysis relevant 
differences between the above two inspiral families are more pronounced for unequal 
mass BH binaries. Incidentally, we also find similar behavior while plotting Xi{q) and 
$(g) for maximally spinning m = 20 Mq BH binaries. In this case, the minimum Xi 
value occurs around q = 3 and this is consistent with the trend observed in figure [2] 
It is reasonable to suspect that the spin-squared additional terms are influential only 
for maximally spinning BH binaries. However, we observe qualitatively similar M{q) 
estimates for binary configurations having moderately spinning BHs (yi = X 2 = 0.75) 
as well as having mildly spinning BHs {xi = X 2 = 0.5). Moreover, the M{q) 
computations indicate that the effect of additional 3PN order terms in d^/dt and 
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dx/dt equations are non-negligible even for single-spin compact binaries. In figure [31 
we plot A$(( 7 ) and Ai{q) while considering I and fc-based single-spin waveforms for 
m = 3 OMq,0i{xo) = 60° binaries in edge-on orientations. We find that these plots 
fairly resemble their double spin counterparts shown in the right panel plots of figure[I] 
Additionally, the variations in the initial values of cos“^(fc • j) and cos“^(Ar • j) are 
similar to the double spin binaries, while varying q value from 1 to 10. Therefore, 
we conclude that it may be beneficial to keep the additional 3PN order terms in the 
differential equations for $ and x while modeling inspiral GWs even from single-spin 
binaries via the precessing convention. 



1 2 3 4 5 6 7 8 9 10 

q 


Figure 2. Plots that display variations in and A4 as functions of q for edge- 
on maximally spinning BH binaries having m = 30 Mq, 40 Mq and 50 Mq while 
keeping 9i{xq) to be 60°. All other initial angular parameters are identical to 
those used in the right panel plots of figure^ The position of the dip plots 

is shifted towards higher q values for higher ra binaries. A possible explanation 
relies on the more influential contributions from the various dot products that 
appear in equation (0 for such binaries having q roughly in the 4 to 9 range. 


The LSC also developed lalsuite SpinTaylorT2 code that employs a slightly 
different version of dx/dt equation while implementing the Tbased precessing 
convention. Note that it is possible to construct various PN approximants to model 
inspiral GWs by using the same PN accurate expressions for the conserved energy, 
far-zone energy flux and the energy balance argument [14] . The differential equation 
for x in the SpinTaylorT2 approximant arises from the following considerations. We 
begin by displaying (symbolically) the 3.5PN accurate expression for dx/dt, employed 
in the lalsuite SpinTaylorT4 code, as 


dx 

dt 


64 c° 
T4 5 Gm 


rjx^ ^1 + Aix + Ai,5 + A 2 x'^ + A2.5 x^^"^ + A^x^ + ^3,5 , 


( 20 ) 


where the coefficients Ai are functions of rj, tt, log(6 x) and 7 e (Euler’s constant) when 
including only the non-spinning contributions. However, the coefficients from Ai s to 
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Figure 3. Plots that compare single-spin SpinTaylorT4 inspiral waveform 
families in the I and fc-based approaches. These A$(q') and A4{g) plots are 
for edge-on m = 30Mq binaries containing maximally spinning dominant BH 
having 9i{xq) = 60°. The qualitative and quantitative nature of these plots are 
essentially identical to the double spin binaries of figure FTI 


A 3,5 additionally depend on Vi, X 2 , (si • Z), (^2 ■ Z), (si • S 2 ), Xi and X 2 to incorporate 
various spin contributions. The explicit expressions for these Ai coefficients may be 
obtained either from equation (3.16) in [28] or from the lalsuite SpinTaylorT4 code 
itself. With these inputs, one defines the differential equation for x in the SpinTaylorT2 
approximant to be 


dx 

dt 


_ 64 

T 2 5 Gm 


64 


5 Gm 


|l + Ai a; + Ai,5 x^G a^x'^ + A2.5 x^G A^x'^ + A3.5 


( 21 ) 


|l + A'^ a; + A'^ 5 x'^G A^x"^ -\- A 2 5 x^G A'^x^ + Ag 5 x'^g'^ 


( 22 ) 


This construction ensures that the coefficients A' are going to depend on various Ai 
coefficients due to the binomial expansion of the denominator of equation (| 2 ip that 
includes all the x'^G contributions. For example, the Ag 5 coefficient is going to depend 
explicitly on all the Ai coefficients whereas the A 2 coefficient depends only on the Ai 
and A 2 coefficients. We observe that the resulting spin-squared terms are usually not 
incorporated into the Ag and Ag 5 coefficients in the lalsuite SpinTaylorT2 code of 
LSC. 

We construct the differential equation for x in our SpinTaylorT2 approximant 
in a similar manner. Therefore, the various A' coefficients are given in terms of Ai 
coefficients which can be extracted from our equations ([9|), (lA.ll) and (IA.2I) for dx/dt. 
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Figure 4. Plots for A4> and Jv{ as functions of q that invoke I and fc-based 
SpinTaylorT2 approximants. We are dealing with edge-on maximally spinning 
BH binaries having m = 30Mq and 9i{xq) = 60°. The A4> and Ai estimates are 
similar to those obtained using SpinTaylorT4 approximant as shown in the right 
panel plots of figure^ 


For example, the expressions for A 2 and A 2 approach are given by 


Ao = 


/34103 , 13661 
U8144 


13661 59 2\ 1 / IN/ I n\ 

‘ 201 ^'^ ~ (247(si • S2) - 721(si • k)[s2 ■ k)j 


+ Xi (2 ~ ~ 

+ ^2 X2 (I (3 (S2 • fe)2 - 1) + T (7 _ (s, . fc)2)) 


(23) 


A'2 — —A2 + A^ 


r/34103 13661 

A18144 2016 


V + - ^» 7 XiX 2 (247(si • S 2 ) - 721(si • k){s2 ■ k)j 


+ xl (3 (si • - 1 ) + F (7 - (si ■ kf)^ 

+ ^2 xl (3 (S2 • A;)2 - 1) + T (7 _ (s2 • fc)^)^ 


+ 


743 1177 


(24) 


In the fc-based SpinTaylorT2 approximant, we incorporate our additional 3PN 
accurate spin-spin terms which contribute only to the Ag coefficient. In other words, 
Ag = — Ag -I- A^ 5 -I- 2Ai A 2 — Af contains our additional 3PN order spin-squared terms 
via the SpinTaylorT4 based Ag coefficient. We would like to emphasize that we are not 
incorporating possible 3PN and 3.5PN order spin-spin terms arising due the binomial 
expansion. This is to make sure that the I and fc-based SpinTaylorT2 approximants 
differ only by our additional 3PN order spin-spin terms. We explore the implications 
of these additional 3PN order terms again with the help of A4{hi, hk) computations, 
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where hi and h^. now stand for the Z-based and fe-based waveform families that employ 
the above described PN accurate expressions for dx/dt. In figure IH we plot Ai{q) and 
A<i>(( 7 ) for edge-on m = 30 Mq maximally spinning BH binaries having dominant 
spin-orbit misalignment equal to 60° at xq. The M and A$ estimates are similar to 
those obtained using SpinTaylorTd approximant as shown in the right panel plots of 
figure [TJ 

It may be argued that the inclusion of our 3PN order spin-squared terms 
to the differential equation for x is not desirable. This is because we do not 
incorporate 3PN order contributions to dx/dt that arise from the next-to-leading 
order spin-spin interactions. These terms are not included as they are yet to be 
computed in the literature. Clearly, the present match computations should be 
repeated when these 3PN order contributions are explicitly available. However, this 
should not prevent one from probing the data analysis implications of the adiabatic 
approximation employed in the X^-based precessing convention of [37] . This could be 
especially useful while invoking PN accurate differential equation for x that incorporate 
higher order spin contributions as pursued in lalsuite SpinTaylorTd and lalsuite 
SpinTaylorT2 codes, developed by the LSC to implement the XN-based precessing 
convention. Hopefully, the data analysis implications of our 3PN order terms should 
provide additional motivation to explicitly obtain the next-to-leading order spin-spin 
contributions to dx/dt. In contrast, it should be noted that the additional PN 
contribution to d^/dt, given by equation ([7]), are fully accurate to 3PN order. The 
fact that such additional precessional effects enter the expression for d^/dt only at the 
3PN order should justify naming our approach as fc-based precessing convention. It 
is worthwhile to remember that precession induced modulations enter the differential 
equation for the orbital phase at I.5PN order in the absence of precessing convention. 

4. Conclusions 

We developed a precessing convention to model inspiral GWs from generic spinning 
compact binaries that employs PN accurate orbital angular momentum X to describe 
binary orbits and to construct the required precessing source frame. The main 
motivation for our approach is the usual practice of using PN accurate precessional 
equation, appropriate for X, to evolve the Newtonian orbital angular momentum 
Xn while constructing inspiral waveforms. We showed that this practice leads to 
higher order PN corrections to = a; equation. A set of differential equations 
and the quadrupolar order GW polarization states in certain frame-less convention 
are developed to model inspiral GWs in our X-based approach. We explained why 
the differential equations for the orbital phase and frequency will have additional 
3PN order terms in our approach compared to the usual XN-based implementation 
of precessing convention. The influence of these additional 3PN order terms were 
explored with the help of match computations involving X and XN-based inspiral 
waveforms for spinning compact binaries with physically equivalent orbital and spin 
configurations at the initial epoch. We adapted both lalsuite SpinTaylorTd and 
LALSUITE SpinTaylorT2 codes, developed by the LSC, while implementing XN-based 
precessing convention. The resulting match estimates indicate that our additional 
3PN order terms should not be neglected for a substantial fraction of unequal mass 
BH binaries. 

It will be useful to pursue our match computations for an extended range of 
the relevant parameter space. The present computations should also be extended by 
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invoking 1.5PN order amplitude corrections to both families of inspiral waveforms in 
the frame-less convention (the lalsuite SpinTaylorT4 and lalsuite SpinTaylorT2 
codes do incorporate such amplitude corrections). This requires us to add certain 
additional 1.5PN order amplitude corrections to the usual expressions for /ix and as 
noted in [44]. Investigating the influence of these additional 3PN (1.5PN) order terms 
in phase (amplitude) while estimating the GW measurement accuracies of compact 
binary parameters will be interesting. It should also be worthwhile to probe the 
influence of these additional terms during the construction of inspiral-merger-ringdown 
waveforms from generic spinning BH binaries with the help of the effective-one-body 
approach. 
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Appendix A. 3.5PN accurate expression for dx/dt 


We list below the regular contributions to dx/dt that we employ while implementing 
our version of the SpinTaylorT4 approximant. These contributions, adapted from |28) . 
read 


dx 

dt 


64 

5 Gm 


TJX < 1 -|- a; 


743 


11?7 


-ka;' 


,3/2 


47 25 

47r - y Sfe - — (Ai - X 2 ) <Tk 


+ x" 


/34103 13661 


V18144 
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2016 


V + - ^^XiX 2 (247(si • S 2 ) - 721(si • k){s2 ■ fc)) 


+ Xi (2 ~ ~ 
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where and are given by 


Sk = Xf xi (si • fe) + -^2 X2 (^2 • k), 
cTk = X 2 X2 (S2 ■ k) - Xi Xi (si ■ k). 


(A.2a) 

(A.2b) 


We would like to stress that in the above equations we have merely replaced I appearing 
in equation (3.16) of |28) by k. This differential equation for x is also employed in 
the usual implementation of the SpinTaylorT4 approximant, as provided by lalsuite 
SpinTaylorT4 code, while using I to describe the binary orbits. 
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